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We perform the numerically exact diagonalization study of the t-J model with nonmagnetic impurities to clarify 
the relation between Zn impurities and the stripes. By examining the hole-hole correlation function for a two-hole 
vl8xVl8 cluster with a single impurity, we find that the impurity has a tendency to stabilize vertical charge 
stripes. This tendency is caused by the gain of the kinetic energy of holes moving along the stripes that are 
formed avoiding the impurity. 
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Nonmagnetic Zn impurities substituted for Cu 
cause dramatic change of the electronic states of 
the high-T c cuprates, since Zn induces spatially 
inhomogeneous spin and charge distributions. On 
the other hand, stripes observed in the cuprates 
have intrinsically inhomogeneous charge distribu- 
tion Therefore, it is interesting to study the 
interplay of Zn impurities and stripes. From ex- 
perimental side, indications of the close relation 
between Zn impurities and stripes have been re- 
ported: The 1/8 hole-concentration anomaly of 
T c , which is an indication of the stripes, is in- 
duced in Bi 2 Sr2CaCu208 and YBa2Cu307_5 by 
Zn substitution [^[jj. 

In order to examine the relation between Zn im- 
purities and stripes, we perform the numerically 
exact diagonalization study of the t-J model with 
nonmagnetic impurities. By examining the hole- 
hole correlation functions for a small cluster with 
a single impurity, we find that the impurity has a 
tendency to stabilize vertical charge stripes. 

We use a two-dimensional t-J model with an 
impurity site. The Hamiltonian reads 
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where c\ a is an electron creation operator at site i 
with spin a with the constraint of no double occu- 
pancy, and the summation runs over all the 
nearest-neighbor pairs. The impurity is located 
at site 0. Here t is the hopping parameter and 
J is the exchange interaction. t lmp (Jj mp ) rep- 
resents the hopping (exchange) strength of four 
bonds around the impurity site. The case of 
timp—t and J lm p=J corresponds to the t-J model, 
while the special case of t imp — J imp =0 describes 
the model with a vacant (nonmagnetic impurity) 
site. We call the latter the Zn model hereafter. 
We note that the paring properties of doped holes 
in the Zn model have been investigated using 
small clusters MM. In the present study, we use 
a \/l8 x%/l8 cluster with a single impurity and 
with two holes under periodic boundary condition 
(PBC) to clarify the relation between Zn impuri- 
ties and stripes. J is taken to be OAt. 

Let us start with the t-J model without im- 
purity. It is controversial whether the t-J model 
itself has the stripe- type ground state ||;@]. F° r 
small-size clusters, the stripe-like state is not sta- 
ble in the ground state but appears as an ex- 
cited state [Q. In the two-hole 18-site t-J clus- 
ter, the ground state is a uniform state with the 
Bi irreducible representation of Ci v point group 
(d x 2_ y 2 symmetry). By examining excited states, 
we identify one of excited states as the stripe- type 
state. This state belongs to the E representation 




Figure 1. (a) V18xV18 cluster with PBC. The 
bold circles represent impurity sites that are re- 
peated by PBC. A part of sites are numbered, 
(b) Hole- hole correlation (n^n^) of the state with 
d X 2-y 2 symmetry in a two-hole 18-site cluster, j 
runs over the sites numbered in (a). Open cir- 
cles and solid squares represent the correlations 
for the t-J model and Zn model, respectively, (c) 
The same as (b) but for the state with p x sym- 
metry. 



having degeneracy of p x and p y symmetries. We 
note that the degeneracy is simply lifted by in- 
troducing anisotropic terms that favor either the 
vertical or horizontal distribution of holes. We 
have confirmed that the excited state identified as 
the stripe state lowers its energy under the pres- 
ence of a potential that stabilizes vertical charge 
distribution and becomes the ground state when 
the potential exceeds a critical value as discussed 
in I- 

To see the charge distribution in the ground 
state and the excited state, hole-hole correlation 
(wjUj') is examined, where (• • •) is the expecta- 



tion value for a given state, n'* is the number 
operator of hole. One of the two holes is fixed on 
the position 7 in Fig. 1(a). The dependence of 
the correlation on the position j of the other hole 
is shown in Fig. 1(b) and(c). In the state with 
d x 2_ y 2 symmetry (the ground state), the posi- 
tions j=X, 3, and 12 show maximum in the cor- 
relation, indicating that the probability that the 
holes are on diagonal sites is greater than nearest- 
neighbor sites p. In contrast, the two holes in 
the state with p x symmetry obtained after lifting 
the degeneracy applying a very small anisotropic 
interaction have large probability to be found as 
a pair on sites with the distance of 2 and y5 (see 
the positions j=0, 4, 9, and 13). It should be no- 
ticed that the correlation along the vertical pair, 
for instance, 7-9 is greater than that of a horizon- 
tal pair of 7-5 under PBC. This means that the 
vertical distribution of the two holes is favorable 
in the p x state. Hereafter we call the p x state the 
"stripe" state, while the d x i_ y i state the "uni- 
form" state. 

The nonmagnetic Zn impurities simulated by 
setting iimp=^imp=0 induce spatially inhomoge- 
neous charge distributions. Figure 2 shows the 
difference of the hole density at site j from the 
average value of hole density (2/18~0.11). It is 
found that the hole density near the Zn impurity 
(site 0) is reduced, while the enhancement of the 
density takes place at sites apart from the impu- 
rity. The loss of the kinetic energy near the impu- 
rity induces such hole distribution ]5|JlO|| . This be- 
havior is seen in both the "uniform" and "stripe" 
states. 

An important effect of the impurity is clearly 
seen when the energies of the two states are 
compared: The energy difference between the 
"stripe" and "uniform" states, -Estripe - -^uniform, 
decreases from 0.34i to 0.15t in the presence 
of the impurity. This means that the "stripe" 
state is more sensitive to the impurity. In other 
words, the Zn impurity has a tendency to sta- 
bilize the vertical charge stripes. More interest- 
ingly, the "stripe" state could be the ground state 
by introducing long-range hopping terms whose 
magnitude is within a realistic parameter region 
jn]]. An intuitive explanation about the sensi- 
tive behavior of the "stripe" state to the impu- 
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Figure 2. Difference of the hole density (AnJ) 
at site j from the value for uniform hole distribu- 
tion (2/18~0.11). The position of Zn impurities 
and site labeling are shown in Fig. 1(a). Open 
squares and solid circles are for the state with 
d x 2_ y 2 symmetry (the "uniform" state) and with 
p x symmetry (the "stripe" state), respectively. 



rity is easily obtained through the examination 
of the hole-hole correlation. From Fig. 1(b), it 
is found that the correlation (ntyrij} in the "uni- 
form" state shows a prominent enhancement only 
at j—3 whose position is the farthest from the 
impurity. Putting this result and other hole-hole 
correlations (not shown here) together, we can 
say that the holes are predominantly confined to 
the region surrounded by the impurities. On the 
other hand, the "stripe" state exhibits enhance- 
ments of (n^n^) at not only j=3 but also j=8 
and 9. Therefore, the holes in this state can get 
ability to gain the kinetic energy through the mo- 
tion in the region between the impurities. This 
contributes to lowering the energy of the "stripe" 
state under the presence of the nonmagnetic im- 
purities. 

In summary, we have perform the numerically 
exact diagonalization study of the t-J model with 
nonmagnetic impurities to clarify the relation be- 
tween Zn impurities and stripes. By examin- 
ing the hole-hole correlation functions for a small 
cluster with a single Zn impurity, we have found 
that the Zn impurity has a tendency to stabilize 



vertical charge stripes. This comes from the gain 
of the kinetic energy trough the motion of holes 
along the stripes that are formed away from the 
impurity. 
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